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mechynism mm elexibly stacked slides 

The mechanism examined is the common slider - crank 
mechanism in which the slider mass is connected with the coupler 
by means of a linear spring. The crank end the coupler are elastic 
and their masses and inertias have also been considered. Equations 
for Kineto-Elasto dynamic inalysis of the Mechanism have been derived 
and an iterative scheme is suggested to solve them, friction is 
assumed small and the effect of clearances in joints and tolerances 
on links have been neglected. Some example problems a re worked out 
and the results are compared with the available ones. 
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CHAPTER I 


1ST ROIC CT ION 


1.1 General Introduction 

luring last fifteen years a lot of -work has teen done in. 
the field of design of mechanisms. Many investigators have contri- 
buted towards kinematic synthesis of mechanisms and it has proved to 
be a powerful design tool for engineers. Almost all the techniques of 
kinematic synthesis of mechanisms* work well and give acceptable re- 
sults for low-speed mechanism. However, in common, they suffer with 
one major shortcoming. This is the "rigidity" assumption which pre- 
vails throughout Hie literature with few notable exceptions. 

In high - speed mechanisms the inertia forces are large 
enough to cause speed - fluctuations as well as appreciable elastic 
deformations in the links of the mechanism, which can not be over- 
looked. The effect of inherent elasticity of a mechanism system may 
not be pronounced in low speed light weight mechanisms, but in high 
speed mechanisms and in particular in high precision mechanisms, the 
analysis may not be perfect unless this effect is considered. Eor 
example a kinematically synthesized gripp er mechanism considering 
all of its links rigid, may miss the target if it is run at a speed 
more than one - half of its design speed. 

In certain machines with reciprocating or oscillating 
motions proper timing of the motion is a prerequisite to smooth 
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operations. At high speeds owing to the large inertia forces., the 
main drive undergoes fluctuation of speed and the link - members of 
the mechanism are elastically deformed. This causes a difference 
between the actual kinematic characteristics of the mechanism at 
certain predetermined positions of the main drive and that predicted 
from the point of view of rigid body assumptions. 

The deviation of the mechanism characteristics from the 
expected ones perhaps motivated the workers to search for some new 
analysis and synthesis techniques in vhich the effect of inherent 
elasticity of mechanism system was taken into account. The analysis 
procedures of the mechanisms have been grouped in the following main 

r H 

categories. 11 j 

Static itoalysis : Determination of Inertia forces, stresses, strains 

and deflections in members and joints due to 
external and / or gravitational loading. 

Kinematic Analysis : Examination of the displacement, velocity ratios, 

acceleration ratios etc., of a mechanism with all 
of its members regarded as rigid. The reference variable is usually 
a position parameter. 

Dynamic Analysis i Determination of the displacements, velocities 

accelerations etc. of a mechanism, including 
derivations of inertia forces of a mechanism made up of rigid 
members. The reference v art able is time. 
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Elastic Analysis : Examination of stresses and- deflections of an 

elastic system due to static load in order to 
determine system flexibilities or stiffness. 

Elastodynanic Analysis : Examination of displacements, velocities, 
accelerations 9 stresses, strains etc., of a moving elastic 
mechanism. Inertia forces are calculated by assuming all the members 
rigid. 

Itineto-elasto dynamic Analysis (KEDA) : Examination of the displace- 
ments, velocities, accelerations, stresses, strains etc., of 
a moving elastic mechanism. Effects of elastic deformation upon the 
inertia forces are included in the analysis. 

An ideal aialysis procedure of high speed mechanisms should 
not only consider the fluctuation# in speed of the main drive owing 
to the inertia forces and the inherent elasticity of the mechanism, 
but it should also take into account the damping characteristics of 
the system, clearances in the joints, tolerances on the links of the 
mechanism, etc. However, such an analysis will be very much complica- 
ted and nobody hitherto has analysed the mechanisms considering all 
these effects together. However, attempts h aye been made to dynami- 
cally analyse and synthesize mechanism with assumptions of varied 
degrees. In some cases particular links were considered to be elas- 
tic and masses to be concentrated at the hinges and so on. The 
analysis and synthesis of slider - crank mechanism seemed to have, 
judging from the published papers cn this topic, attracted least 
attention of all the planar mechanisms. 
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1.2 literature Survey 

Ihring past few years valuable contributions have been made 
in the field of dynamics of mechanisms. Perhaps, due to the inaccurate 
response of the kinematically synthesized mechanisms at high speeds, 
the literature has recognized the need for dynamic analysis and syn- 
thesis techniques. However, in a great majority of work the dynamic 
analysis has been performed considering rigid links. 

Early work in dynamics of mechanisms concerned themselves 
with deriving the velocities, accelerations etc. of mechanisms. 

Quinn | 31, 32~j has pioneered in the energy distribution method and 
has set the stage for further work. Eor determining dynamic charac- 
teristics of mechanisms, Hiirchhom j~l6j has given the ’’rate of change 
of energy approach” in which the time rate of change of total energy 
of the median! cal system is equated to the power input. Thus, if the 
power input function is known, this method will yield the time res- 
ponse of the system. 

Elastic deflections of significant magnitude in high speed 
mechanisms due to their inherent elasticity have attracted the wor- 
kers and various techniques ‘have been developed to account for this 
effect. The use of springs in the mechanisms h we also been made in 
order to alter the dynamic performance of the given mechanism in 
high speed motion. This may improve the response of the system, 
reduce the shaking forces transmitted to the fraae or support or 
produce a performance vhich would otherwise be unsatisfactory with 
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have dealt with "elastic 


the rigid member mechanism. Van Sickle and Goodman | 4-3 j used springs 

to adjust the speed of mechanisms. Imam |^l9*j examined the response 
of a four - bar function generator with a longitudinal spring as a 
part of its coupler. Davidson 5 j has examined a slider - cranio 
mechanism with the slider mass connected by a spring. In J~40~j Benedict 
states that "the addition of springs is one of the easiest ways to 
balance torque and obtain a smooth - running machine". 

Some aithorsj^2, 3, 19, 20, 25 

complex" systems i.e. a mixed mechanism in which few members are 
elastic and the remaining are taken to be rigid. Lagrsngian mechar- 
nics or Energy methods are used to derive the equation of motion. 
Generally, as the solution of these problems is tedious, only one 
member is taken to be elastic and that too with only one elastic deg- 
ree of freedom. 

Lubowsky and Freudensteinj^6 J have made notable contribu- 
tion to the study of backlash aid have modeled an elastic mechanical 
joint, which they have called as "impact pair", and have determined 
the dynamic equations for the relative motion of such joints. 


Slider - crank mechanism has been the centre ox study of 
several workers. Jasinski, Lee and Sandorp20, 21 j have studied the 
stability of this mechanism having a flexible connecting rod* 

Henbane r, Cohen and H all £ 2 9~] have made a notable contribution by 
examining the vibrations of the elastic connecting rod of the slider — 


crank mechanism 
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liniecki 26 | considers the fluctuation of input crank 

— i. 

speed due to the effect of inertia forces. He deals with a procedure 
for the kinematic redesign of the slider - crank mechanism so as to 
miniimi.se the error between the actual and the desired slider veloci- 
ties at three predetermined positions of the crank. 


Davidson 


CO 


has examined a slider - crank mechanism connec- 
ted in series with a spring and mass. He has taken crank and coupler 
to be massless and rigid and the only inertia load is due to the mass 
of the slider. He has also assumed a constant driving speed and with 
these assumptions he has synthesized the mechanism as a function 
generator. 

Sheerwood J~30 considers the dynamic synthesis of a slider- 
crank mechanism so as to minimise the error between the actual and 
the desired time taken by the slider to travel a prespecified dis- 
tance. Sheerwood aid Hockey |j34 , 3Tj have considered the redistri- 
bution of coupler mass of a slider - crank mechanism from balancing 
point of view. 


The literature reviewed above has not provided a method for 
"Kineto - Elastodynsmic .Analysis" (KEDA) of mechanisms in viiich all 
effects are considered. Only recently some contributions hare been 
made towards KEDA of mechanisms. 


Erdman and SandorQ7~j have called KEDA to be "a frontier" 
in mechanism design. Sadler and Sandor j~ 33~j used harmonic analysis 
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to analyse the kineto - elastodynamic deflections of the extension 
of the coupler of a four - bar path generator linkage. Erdman f 
Imam, Sandor and Oakberg J~8, 9, 10 s lijhgve dealt with a general 
method for KEDA and Kineto - Elastodynamic Synthesis (iCEDS; of 
mechanisms using Kineto - Elastodynamic Stretch Eotation Operator 
(KEDSRO). They have started v;ith the kinematic analysis to get 
system forces and elastic deformations 5 derive input link accolera r- 
tion using law of conservation of energy and get velocities and 
accelerations of remaining links from what they call "Kineto - Elasto- 
dynamic family tree", obtained by differentiation of KEDSRO. The 
entire procedure is iterated till a suitable convergence criterion 
is satisfied. Imam [~ 18j deals with a KED.6 and KEDS technique and 
without using "KEDSRO" synthesizes a mechanism with minimum weight 
criterion. 


1.3 Aim Of The Present Work 

As stated earlier kineto - elastodynamic analysis of high 
speed mechanisms is a recently developed field of study and deserves 
attention of mechanism designers. Particularly in high speed high 
precision mechanisms (e.g. gripper mechanisms of computers and prin- 
ting machines etc.) the elastic deformations are large enough to 
affect the mechanism response and hence the desired accuracy can 
only be achieved if the effect of the inherent elasticity of the 
mechanism is also considered during analysis. 
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In the present work kineto - elastodynamic analysis of the 
eccentric slider - crank mechanism with flexibly connected slider 
has been performed. All the links have been considered to be elas- 
tic. The masses and inertia effects of all the links have also been 
considered. However, the effect of the clearances and tolerances, 
damping characteristics of the mechanism and the friction effects 
at the junction of two links have been neglected. 

The common slider - crank mechanism is a particular case 
of the mechanism considered here when the stiffness of the spring 
connecting the slider is considered to be very hi^a (theoretically 
infinite ) . 




CHAPTER II 


FORMULATION OP THE PROBLEM M D SOLUTION TECHNIQUE 

2.1 Formulation Of The Problem 

Figure 1 represents the mechanism under study. The stiff- 
ness k of the spring connecting the slider with the coupler is 
assumed to be linear. The crank and coupler are assumed to be 
linearly elastic and their elastic and section modulii are taken to 
be uniform and constant. The mass of the coupler is concentrated at 
its centre of gravity G (which is taken at its mid - way for 
simplicity.) The gravitational forces, being much. small in compari- 
sion to inertia forces, are neglected. The joints are assumed to be 
frictionless and slider - friction has also not been accounted. 

There is no play in joints. Linear superposition of rigid - body and 
elastic deformations is assumed. 

The mechanism driver rotates at an average angular velocity 

Jfl. 'The main drive undergoes fluctuations of its speed due to 

inertia of crank, coupler and the slider mass and at the steady - 

state the input crank angular velocity © is a function of input 

angle © only. The distance X represents the actual location of 

c 

1 

coupler end point C while X represents location of point D. 'When 
the mechanism is motionless point D on the slider coincides with 
point G on the coupler. Por very large values of k (a rigid connec- 
tion between slider aid coupler) C and D would coincide for all 
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phases of the motion. But in the mechanism under consideration the 
slider is flexibly attached and has mass n^. Therefore point D is 
displaced by an amount (x - JC ) with respect to the point G. The 
displacement 2 ' would be the algebraic sum. of "rigid" displacement 
and elastic deformation of the coupler end C in x - direction. 

2.1.1 The Analysis Procedure : 

Since crank and coupler are having mass and inertia and 
are elastic, it is not possible to obtain displacements, velocities 
and accelerations of different links at various positions of the 
crank directly. The displacements would constitute of two compo- 
nents - the dynamic displacements and the elastic deflections. The 
elastic deflections cannot be obtained directly as they are func- 
tions of system forces. Ttor getting 'true' elastic deflections, 
’true’ system forces should be known. But the ’true’ Inertia forces 
depend on 'true' accelerations which In turn depend upon 'true' 
elastic deflections. Thus, system forces and elastic deflections 
are inter - dependent and hence a suitable iterative scheme is re- 
quired to evaluate elastic deflections. The following procedure of 
analysis is adopted. 

The equation of motion of the slider and the energy ba- 
lance equation (the equation obtained by equating rate of dhange 
of total energy to the power input to the system) constitute the 
equations of analysis. The ordinary differential equations thus 
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obtained are non - linear aid coupled. The values of velocities and 
accelerations of centre of gravity of the coupler and displacement 
of tiie coupler - end G (Fig. l) are required in order to solve them 
to yield the velocities and accelerations of the crank and slider, 
and the slider displacement. Since the required quantities them- 
selves are functions of velocity and acceleration of the crank, 
these equations cannot be solved directly. An iterative scheme is 
adopted to solve them. 

Firstly, all the velocity aid acceleration terms required 
to solve the above two equations are kinematically expressed in terms 
of velocity and acceleration of the crank. The displacement of the 
coupler - end C is aLso expressed kinematically. The equations are 
then solved together using fourth order Range - Kutta method, yield- 
ing dynamic velocities aid accelerations of the crank and the slider, 
aid the slider displacement. The system forces are then obtained. 
Elastic analysis via flexibility approach yields element deflec- 
tions. Then the velocities aid accelerations of elastic deforma- 
tions may be obtained. Knowing accelerations of deformation the 
dynamic accelerations aremodified. The modified system forces, ele- 
ment deflections etc. may now be obtained and the process is iterated. 
This scheme is iterated to desired accuracy so that system forces, 
elastic deformations aid crank accelerations are best modified. 

Sow the modified values of velocities aid accelerations 
of centre of gravity of the coupler and coupler - end displacement 
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may be obtained. With these modified values the differential equa- 
tions of analysis are again solved to give a better approximation 
for the velocities and accelerations of the crank and the slider and 
the slider displacement. The entire procedure is now iterated till 
the mechanism characteristics obtained from two consecutive itera- 
tions match to the desired accuracy. 

2.1.2 Derivation of Equations of Analysis 

Refering to fig. 1 the differential equation of motion of 
the slider may be written as 


• * 

m X + k (X - Z ) = 0 

1 c 

Or Z + w 2 Z = w 2 Z (2.1) 

n n c 

where w = . 

n ym 


from kinematic considerations, one gets 



Gos 

0 = 

Hi - /-a- 

! ! i 

' 4* 

r 

1 

■) ^ 

Sin 8 t 

1/2 




L_ t 

! 



^ *"*“"* J 

■v 2~ 

and 

X 

= r Cos 8+1 

i 

- j 

-f- Sin 8 + £ 

i 

c 



1 

1 1 

j 




where 

r = radius of crank 

1 = length of connecting rod 

a = amount of offset 

8 = angle of the crank with x - axis 
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0 = angle between the coupler and x - axis 

X = "rigid" or kinematic displacement of coupler end C. 
o 

It has been shown. jj26j that if the sum of y and -r 
is snail (e.g. less than 0.3) 


1 - ( ~ Sin 0 



with an error of nearly 0.08%. 


I 1 - \ ( f Sto 6 + f ) 2 

L - 


can be re writ 

xen as 

= r Cos 6 + 1 

1 ( p a, lel 

1 - i 1 i Sin e f f f 




( 2 . 2 ) 


The working section of the torque - speed characteristics 
of a three phase a/c motor may be successfully approximated by a 

parabolic, equation when referred to the main shaft f~ 26~| . 

.2 

. — » 

1 


a„ - & 2 0 


(2.3) 


where 

T 

& a 2 

§ 


external torque overcoming fluctuations in speed 
referred to the main shaft. 

Coefficients of applied torque 
Crank velocity 


Applying the Energy Balance to the mechanism system yields 

d ( T.. E.) = H0 (2.4) 

Where, T.E. is the total energy of the system at any phase and 
includes the kinetic energy of crank, coupler and slider mass, the 
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potential energy stored in the spring and the Strain Energy stored 
in crank and coupler owing to elastic deformations. But due to the 
presence of a spring, the strain energy stored in the crank and the 
coupler will bo much smaller in comparision with the kinetic and 
potential energies and hence is not accounted. 

Therefore, 

T.- E. = IC. E. + P. E. 

where, 

K. E. = kinetic energy of the mechanism 
P. E. = Potential energy stored in the spring 

Substituting in equation (2.4) results 

d (K.E. + P.E.) = T a e 


Or 


-A 

d0 


(K.E.) + (P.E.) = T 


2 

» 

a 1 “ 6 


(2.5) 


The K.E. of the mechanism is given by 

2 2 *2 • ^ 2 

1 _ * _ v + eu + m 0 Y„ + m A 

K.E. = -^ I © + I p $ 2 & 2G- 1 

2 0 «■ 


( 2 . 6 ) 


vtfxere 


I 

o 




0 


X & and Y & 


Moment of inertia of crank about the satis of 
rotation 

Moment of inertia of coupler about its centre 
of gravity 

Angular velocity of rotation of the coupler 
x and y - directional velocities of centre of 


gravity of the coupler 





17 

X = linear velocity of tie slider 
= Mass of the slider 
nig = Mass of the coupler 


Differentiating equation (2.6) with respect to 8 and setting 

* 2 

Y = © , one gets 


i ( K.E.) - |l 0 T , +? f 


• • • 


I Q 0 0 + n 2 (X^ X G + Y & Y fl ) 


.rx' 2 / 

* i 


I !t‘ 

+ 2Y X X 


(2.7) 


Since 


■i (x 2 ) - x ’ 2 § ♦ ax x’ x’ 


aid 


■£ ( ?) - H* 


so on. 


© 


The primed quantities in the above equations represent 
derivatives with respect to © . 

The potential energy stored in the spring is given by 

P.E. = l \ ( X - X ) 2 
h. G 


( 2 . 8 ) 


Differentiating (2.8) with respect to © ? one gets 

” (P.E.) = k (X - X c ) (X* - X G ) 

Substituting (2.7) end (2.9) in (2.5) and simplifying, one gets 

'2 


(2.9) 


r- I m X 


i tr 


^ - Y (ag + m 1 X X ) - 


I G 0 0 + m 2 (X & X & + Y ft Y g ) 
- k ( X - x c ) (x' - t ' c ) 


( 2 . 10 ) 
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Since 

*• 

Z 


*2 » ** t 

e z + © z , 


equation (2.l) yields 

X ‘’ = Y C \ (* c - X) - ^2- ) (2.-11) 

Pnom (2.1l) aid (2.10), one cai get 


2 _ 

I 


- & 2 Y - 


_L 

v/T 


l * •• 

1 ^ 0 0 

L G 


+ m 2 \ + \ V ; 


- k (Z - Z G ) (z' - Z G ) + m 1 z' (Z - Z c ) 


( 2 . 12 ) 


Equations (2 . 1 1 ) and (2.12) constitute the equations of analysis aid 

are two non - linear coupled ordinary differential equations whose 

• « * • 

solution would give the values of X, Z, ©, Z and © . 


2.1.3 Elastic Deformation Of The Mechaiism 

Any mechanism may be considered a structure if its rigid - 
body - kinematic degrees of freedom are removed. At every phase the 
mechanism is "frozen" by removing its rigid - body - kinematic degrees 
of freedom and then its analysis as a structure yields the elastic 
deformations. 

A mechanism is composed of various combinations of elements 
each of ’Aida can be represented by a known structural model. The 
deflections of the entire mechanism system may be derived by per- 
forming an elastic analysis using flexibility approach. The mecha- 
nism will have system, or generalized external forces acting upon it 
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P.! 3 -1, 


, m, 


which will be represented by the column matrix 
where m is the number of system forces. The number of the elastic 
degrees of freedom n of the system is the sum of elastic degrees 
of freedom of its elements, each degree of freedom being represented 
by an element coordinate. For the mechanism under consideration the 
number of elastic degrees of freedom n is taken to be 7 while 5 
generalized forces act on the mechanism. (Fig. 2 and Fig. 3) 


2. 1.3.1 Matrix Of Element Flexibilities 


The mechanism constituting of crack m d coupler only is to 
be analyzed for elastic deformation at each phase. It has got only 
one rigid - body - kinematic degree of freedom. Therefore, it may 
be converted into a structure by modeling the Input link as a canti- 
lever or fixed - free beam. As the end C of the coupler is assumed 
to move along x - axis only, the coupler is modeled as a bean hinged 
at the end A and having double roller support at the end C as shown 


in Figure 2. 


The element flexibility matrix. F cri of the crank is given 


by j^8 ] 


L 


r ' 

1 F 

t cr 


r/A o E 1 
0 
0 


r/3E 1 I 1 
r/2E 1 I 1 


r/2E 1 I 1 
r^i I 1 


(2. 


where 


A = Area of cross - section of the crank 

0 

E = Modulus of elasticity of the material of crank 

1 

1^ = Section - modulus of the crank. 


13) 
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The element flexibility matrix of the coupler can be cons- 
tructed by using C astigliano ’ s Theorem. 

Re fe ring to PIG. 4 the strain energy per unit volume stored 
in the material of coupler is given by 


. 1/2 2 

1 (fj_ Cos 0 + R c Sin 0) 


i) 1 (f 4 Cos 0 + R c Sin 0 - f g Oos 0 + f ? Sin 0 )‘ 


^ 2 (P* Cos 0 - f. Sin 0 s2 


E 2 h 


2 dz 


f^ z Sin 0 + z Cos 0 + f g (z ~ -|) Sin 0 + - 


E 2 I 2 


~)cps-.-0 - f 5 ) 


(2.14) 


UShere 


R , the 
0 


Aq = Area of cross - section of the coupler 
E 2 = Modulus of elasticity of the coupler 
I = Section - modulus of the coupler 
R = Reaction at C 

reaction at the point G, can be evaluated by taking moment of 


all the element forces about A. Then one gets 


fg f 5 £y 

— tan 0 + i Oos K + f 4 t - — 


(2.15) 


I 
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How applying Castigliano 1 s Theorem, the differentiation of equation 


(2.14) with respect to f results 


r bu 

O, = r-r- = 




r^ 2 (f^ Cos 0+E c Sin 0) dz 

\ -m a 


I 

u 0 


E„ ^ 


Cos 0 


A 
1/2 


^ (f , Cos jZ? + E Sin 0 - f ' Cos 0 + f ' Sin 0) dz 


E 2 Aj 


Cos 0 


which results after simplification as 


1 Sec 2 0 tan 0 Sec / 1 Sec 2 / 

O = f + f f 

4 4 E 2 ^ + J 5 E 2 6 2 E 2 ^ 


(2.16) 


Differentiation of equation (2.14) with respect to f g results 


r 

6 c 


^ U 

~b~z* 


r/ 2 (f^ Cos / + Eq Sin /) tan / 


E, 


2 A[J 


dz 


C (f 4 Cos 0 + E q Sin 0 - f g Cos / + f ? Sin /) tan / 


\J 


1/2 


^2 


dz 


fC^ 2 (E Cos 0 - f Sin 0) 2 

i C 4 z 

+ v ipq T _ 32 


' 1 (-f . s Sin {i t IL ! Cos 0 t 1, (z - -i) Sin 0 + f_ (z - hcos ft - f.) 


1/2 


e 2 i 2 


( f - 1 ) dz 


which, on simplification, can be written as 
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C _ f .Sec 0 tan 0 

5 4 E 2 5 


1 tan 2 0 

-f — — 1 


12 B 2 I 2 l E 2 


Sec 0 tan 0 
2 


6 2 E„ Ij, 


(2.17) 


Differentiating equation (2.14) with respect to f r ? one gets 


R = ^ 
6 2^6 


f^ 2 (f^ Cos 0 + Sin. 0) Sin 0 t ax 0 


0 

1 


E 2 \ 


dz 


(f^_ Cos 0 + Sin 0 - fg Cos 0 + Sin 0) 


1/2 


5 2 ^ 


(Cos 0 + ms J^ss J . ) d2 

(- z Sin jZf + Bp z Cos 0 + fg (z - — ) Sin 0 + f ^ (z ~ -|)Gos 0 - f^) 


1/2 


E 2 I 2 


Sin 0 dz - ^ 

J . 


V 2 (IL Cos 0 - f Sin 0) Sin 0 


e 2 i 2 


-g — z dz , 


Or 


S 6 - 


1 See 2 i . Sec 0 tan t . „ 1 (Sec 2 $ + Cos 2 20 

Fi^- f 5 2.^ + f 6[_ 4B 2 ^ 


3 _ . 2 ^ 

. 1 §in_0_ . f 

48 E 2 I | 7 


X Sin 2 0 1 tan 0 (l + Cos 2 0) 

96 E 2 I 2 8 e 2 ^ 


(2.18) 


Differentiating equation (2.14) with respect to f^, one gets 
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S r 


~&f n 


TJ 


1/2 


(f 4 Cos 0 + Bq Sin 0) Sin 0 


E.. 




dz 


0 (f 4 Cos 0 + E Sin 0 ~ f Cos 0 + f 7 Sin 0) Sin 0 


* J 

1/2 

1/2 


E 2 Aj 


dz 


Cos 0 - f Sin 0) 

z Cos 0 dz 


2 E 2 h 


+ l 

1/2 


(- f 4 z sin 0 + ^ z Cos 0 + f g (z - -|) Sin 0 + f ? (z - -|)Cos 0 - f,.) 

— 


I?--: A) 


Cos 0 dz 


Simplification of the above equation results in 


& - 


l 3 Sin 2 0 1 Sin 2 0 

96 E 2 *2 8 E 2 V 


+ f r 1 sin2 ^ * i 5 cos 2 0 


7 i 4 E 2 % 


48 E 2 X 2 


(2.19) 

The element deformations are represented by a n - dimensional 
column matrix and are related to the element forces j~ f^J by the 

relation 

~S = [ ? e ~! ~f ( 2 . 20 ) 

where is n - dimensional column matrix representing the 

element forces and J~ fJJ is the element flexibility matrix (n x n) 
of the mechanism and can be constructed by using equations (2. 13), 
(2.16), (2.17), (2.18) and (2.19) as 
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°1 

0 

0 

0 

0 

0 

0 ! 

0 

° 2 

°3 

0 

0 

0 

! 

0 

I 

0 

°3 

°4 

0 

0 

0 

i 

0 

0 

0 

0 

°5 

C 6 

C 7 

0 ! 

0 

0 

0 

C 6 

C 8 

C 9 

0 

1 

0 

0 

0 

°7 

c 9 

° 1 C 

°11 
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0 

0 

0 

0 

C 11 

C 12 


( 2 . 21 ) 


where 


A o E 1 


C„ = 


r 3 


3E 1 X 1 


CL = 


2 E . ~L. 
1 1 


E 1 E 1 


CL = 


1 Sec fS 

E„ JL 
2 U 


_ ten 0 Sec 0 
6 = E 2 ^ 

r _ iSecLi- 

7 “ 2 E 2 ^ 

1 tan 2 0 

8 " 12 E 2 I 2 1E 2 ^j 

Sec 0 tan $ 

°9 "”2 E 2 Ag 
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10 


11 


12 


1 (Sec 2 $ + Cos 2 fiO l 3 Sin 2 0 
4 E 2 48 E 2 I 2 


u Sin 2 g f 
96E 2 I 2 


1 Sin 2 jZf 

8E 2^ 


1 Sir, 2 2? . I 3 Cos 2 0 

4 E 2 ^ " 48 E 2 I 2 


( 2 . 22 ) 


2. 1.5.2 Eorce Transfer Matrix 

In order to transfer the system forces J P.l j = 1, . m 

; 3 • 


f^i i = 1 , . . . . , n an n x m force 


into element or internal forces 
transfer matrix j Dj is derived by the methods of static analyst 

The transformation relation is given by 


sis. 


— .>■ 
f 


- n "i 


D 

L J 


(2.23) 


The force transfer matrix Dj also relates the system deformations 
represented by m - dimensional column matrix j d . j and element 

, O • 

deformation S K . i end the relation is given by 

"7 = [d] T (2.24) 

Equations (2.20), (2.23) aid (2.24) can be used to get the relation- 




ship between system deformations j d. t and system forces j P. 

Li! L i 

Thus substitution of (2.20 ) and (2.23) in (2.24) results in 


or 


d 

*'*“J5* 

d 


i T 

r P i 

i 

-t 

1 e 


D j P 




(2.25) 

(2.26) 


vile re 
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is a m x m matrix and is known as system flexibility matrix. 

The force transfer matrix jjDj can be constructed by using 
the methods of static analysis. 


Refering to PIGUEES 2, 3 and 4 representing force diagrams 
of the elements, one gets 


f = P 
3 1 


f. = *A 

4 4 


f. 


(2 k .'28) 


H Cos © - R Sin © 
a a 

H Sin © + R Cos © 
a a 


(2.29) 


The system forces cm be written in terms of accelerations 


as 


= I® 

1 o 


1 * 
7 i y 
2 o 


- 

• » 

= 3L 


2 G 

h (X - X Q ) 


(2.30) 


m 2 Y G 


The reactions H and R cm be obtained using free- body 
a a 


diagram of the coupler PIG. 4* Thus, 


a 4 o 
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end R = 
a 


x 5 x 7 

ter. 0 + -f Cos 0 + f ten 0 + ~~ 


Substitutimg (2.31J in (2,29) end simplifying, one gets 


f l 


Sin 9 


'2 1 Cos 0 ' "3 

"4 


+ J ( Sto 8 tm 0 _ Oos e) 


-** |— 

+ P^ (Cos 9 - tai 0 Sin — Sin 9 


(2,31) 


and 

f 3 


Cos 9 


'2 1 Goo 0 ^3 (Sin 6 + 


Cos 9 tan 


4 


±V 

+ P^_ (Sin 9 + tan 0 Cos 9) + — Gos 9 
from (2.28) and (2.32), one can obtain that 


! d ! = 


0 

0 

1 

0 

0 

0 

0 


D 1 

0 

0 

1 

0 

0 


D 6 

0 

0 

0 

1 

0 


3X, 


0 

1 

0 

0 

0 




D, 


8 


0 

0 

0 

0 

1 


(2.32) 


(2.33) 


D 


IX- 


D, 


- Sin 9 


D, = - 


1 cos 0 

Sln . e^tgLi _ Gog e 

Cos 9 - ten 0 Sin 9 
Sin 9 


' 2 


where, 
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Dr- = 


D r 


D, 


7 


3V 


Cos 9 
1 Cos 0 

_ Cos 9 tai $ \ 
- (sm e + & ; 


Sin 9 + ten 0 Cos 9 
Cos 9 


(2.34) 


“B 2 

2. 1.3.3 System Flexibility Matrix 

Finally the system flexibility matrix, an (m x m) matrix, 
may be derived by substituting equations (2.2l) end (2.33) in (2.27). 


where 


f, = 0, 


h = S 1 ? 


i. 


°3 B 6 


*4 = °3 E 7 

S 5 = °3 4s 

*6 - c i d i +0 2 I) 5 + 0 


8 


F 7 = C 1 D n D 2 + C 2 D 5 D 6 + C 9 




*1 

F 2 

P 3 

F 

4 

"5 



? 2 

*6 

F? 

F 

8 

“9 

F 

s 

: “ 

P 3 

*7 

F 

10 

F 

11 

F 

12 

_J | 








*4 

*8 

F_ 

11 

F 

13 

F 

14 



F 5 

"9 

F 

12 

F 

"14 

F 

15 


(2.35) 
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F 

8 


°i 

B 1 

B 3 + 

°2 

D 5 D 7 + 

*9 

5= 

°1 

B 1 

B 4- 

°2 

D 5 D 8 

P 10 

= 


4 

+ °2 

d 2 

0 

+ °10 

'll 

=5 

°i 

h 

B 3 + 

=2 

D 6 B 7 + 

(M 

v— 



B 2 

B 4 + 

°2 

D 6 D 8 + 

F 

= 

G 

Dll 

+ c 

D 2 

+ C 

13 

* 


1 

3 

2 

7 

5 

P 

s- 

c 

D 

D + 

G 

D„ IL 

14 


1 

3 

4 

2 

7 8 

F„ c 

= 

c 

D? 

+ cl 

D 2 

+ G 

15 


1 

4 

2 

8 

12 


(2.36) 


2. 1.3*4 System deformations and their velocities and accelerations 


Having derived the system flexibility matrix, the system 


or generalized deformations 
directly from equation (2.26) 


d. 

3. 


3 = 1, , m may be obtained 


Ref c ring to FIG. 3 one can ccnclade that 


0 = d„ 

e 1 




X G - d 4 
e 


G 


(2.37) 


whore suffix c denotes the clastic deformations (e.g. X n is the 

e 

elastic deformation at point G of the coupler in x - direction) 1 .' 




30 


The velocities and accelerations of deformation nay -then 
be obtained by using finite difference technique. Thus, one nay 
derive 


• t • 

e = e & 

G G 

• t 

0=0 9 

e e 




i • 


Y G 6 
e 



(2.38) 


and 


* * ! • • !! *n 

0 = 90 + 00 =: 

e e e 

. . t •* it 

0 ~ 0 0 + 0 Y 


t «* tt 

9 0 + 9 Y 

e e 


.. , .. „ 

X G = X B 9 + X S T 

e e e 

.. t •* it 

Y s = y g 9 + y g y 

e e e 


.. , .. „ 

X g = Y G 0 + Z c Y 
o e e 


(2.39) 


2.2 Solution Technique 


After having mathematically formulated the problem, the 
next task is to devise some iterative technique to solve the equa- 
tions. As discussed earlier the equations of analysis cannot be 
solved directly. 


Rcfering to equations (2.1l) and (2.12), one may con- 
clude that this set of two non - linear coupled ordinary differential 
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equations may be solved directly to yield the values of X, X, ©, 

*• *• 

X, © etc. provided the displacement, velocity and acceleration 

* • • 

terms (o.g. X & , X^otc.) appearing on the right hand side of 

the equations are known. Hit these displacement, ■velocity aid 

acceleration terms are not only kinematically related to the un~ 

• • • 

known input crank velocity 0 and acceleration ©, but also depend 
upon the elastic deformations, velocities and accelerations of 
deformation. Elastic deformations, velocities and the accelera- 
tions of deformation depend upon the inertia forces which in turn 
depend upon the acceleration terms. Since the acceleration terms 
and tie inertia forces are rater-dependent, this set of equations 
can be best solved iteratively with the set of equations (2.26), 
(2,35), (2.37), (2.38) aid (2.39). 

2.2.1 Iterative Schemes 

As stated earlier the equations of analysis alcngwith 
related equations cannot be solved directly end some iterative 
scheme is to bo adopted for their solution. The procedure adopted 
for their,' solution is shown on the flow chart (Appendix a). The 
steps arc os follows s 

STEP 1 s The "rigid" displacement, velocity and acceleration 
terms appearing cn the right hard side of the 
equations (2.1 1) md (2.12) are expressed in terms 
of velocity and acceleration of the crank. These 
quantities are then substituted in equations of 
analysis. T© start with, the links are considered 


t© be rigid. 
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STEP 2 


. Equations (2.1l) ana (2. 12) can now be solved using 
some numerical technique (e.g, fourth order Bunge - 
Kutt a Method) . Since the above equations constitute 
a boundary - value problem, following procedure is 
adopted for their solution. 

(i) Following values are assumed in the first step. 

(a) Initial guess values for slider displacement 

1 

X, its derivative X with respect to inde- 
pendent parameter © and initial velocity of 

* *2 

the main drive © (or for Y = © ). 

(b) Applied torque constants end a^. 

(ii) Using Runge - Kutt a method values of X, X and Y 
are computed for each 2° of input - crank rotation. 

(iii) The initial values of X^ X and Y at © = 0 and 
their values at © = 2 J\ are compared. In the 
steady state they should be equal, for this de- 
notes that the mechanism, maintains an average 
speed in operation. If X (© = 0) ^ X (© = 2 7Y ) 
the initial guess forX at © = 0 is corrected. 
Similarly if X (© = 0) ^ X (© = 2 Tf ) the ini- 

1 

tial guess value of X is corrected. If 
Y (© = 0) / Y (© = 2 71 ) then the a. term of 
the torque expression (T = © ) is sui- 


tably charged. 
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(iv) The average computed value of Y is compared 

n 2 C' 

with where 3 X. is the desired average 

angular velocity of crank. The averrge value of 

Y is obtained as 31 arithmetic mean of the 180 

computed values of Y based on equal step in ©. 

2 

If average vohra of Y is not close to —i . the 
initial guess value of Y is modified. 

(v) If the requirements of steps (iii) end (iv) are 
not satisfied, with the modified values of X, 

1 , Y and a„ at 8 = 0 the equations (2.1 1) and 
(2.12) ore again solved. The iterative scheme 
is continued till requirements of steps (iii) 
and (iv) are met. 

STEP 3 : Having solved equations (2.1l) end (2. 12) the initial 

position of the mechanism (© = 0) is set. 

STEP 4 : Solution of step (2) yields slider displacement, 

velocities and accelerations of the slider aid the 
crank. Knowing accelerations the system forces are 
evaluated using expressions (2.30). 


STEP 5 : Elastic deformations may, then, be computed using 

expressions (2.26). 

STEP 6 J Using expressions (2.38 end 2.39) velocities end 
accelerations of deformation are computed. 
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STEP 7 


STEP 8 


The values of acceleration terns are modified using 
results of step (6). 

Values of new and old accelerations are compared. If 
the values do not match to the desired accuracy the 
entire procedure from step (4) and onwards is repeated. 
Finally the modified values of inertia forces and acce- 
lerations are obtained. 


STEP 9 : All the velocity terms are modified using results of 

step (6). Displacement of coupler end C is also modi- 
fied using results of step (5). 

STEP 10 : The process from step (4) to step (9) is repeated for 
all mechanism positions between 0 aad 2 JT with a step 
increment of 2°. 

\ 

STEP 11 : If the procedure from step (2) to step (lo) is iterated 
once only, step (2) is ^ain undertaken with modified 
displacements, velocities and accelerations aad the 
procedure is repeated. 


STEP 12 j lew and old values of slider displacement X are com- 
pared. If the values match to the desired accuracy 
the process is terminated. Otherwise the entire pro- 
cedure from step (2) to step (ll) is iterated. 





CHAPTER III 


KESLTIES AND DISCUSSIONS 


To illustrate the application of the mathematical model 
and solution technique developed in Chapter II three numerical 
examples have been vroifced out. All the three examples were taken 
from the piblished literature so that the results obtained may be 
compared with the existing ones. In the first example problem 
the mechanism considered by Davidson 5 J was studied. In other 
two example problems the mechanisms synthesized by Idniecki j 26 
were worked out. 


3.1 Example 1 : 

The mechanism studied in this example problem is sche- 
matically identical with that examined by Davidson [ 5 j . In 
Davidson's model the radius of crank, length of connecting rod 
and the amount of offset are taken to be 2.68 in., 9.38 in. and 
0.0 in. respectively. The driving crank rotates at an average 
speed of 200 rev/min aid the slider weight is taken as 4 lb. wt.. 
The stiffness of the spring is taken to be 42.5 lb/ict. . The 
masses and elasticity of the crank and the connecting rod are 
neglected in his model. Whereas in the present dissertation they 
are taken into account. The crank is assumed to have an area of 

cross - section of 1.0 sq. in. aid moment of inertia about its 

2 

axis of rotation as 1.0 lb-in-sec . The coupler weighs 2.0 lb. 

2 

and has an area of cross - section as 0.75 in. . Its moment 
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Fig. 6 Variation of slider displacement and Rigid 

displacement of coupler-end C with respect 
to crank angle 
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Fig.n Variation of elastic deformation at coupler 
end C with respect to crank angle 
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of inertia about the transverse axis passing through its centre of 
gravity is taken as 0.153 lb-in-sec. 2 . The modulus of elasticity of 
the material of crank and coupler is taken as 30 x IQ 6 lb/in 2 . The 
constant 82 of the motor torque is taken as 1.8. ! 26 | 

3.1.1 Results and Discussion for Example Problem 1 

With the above input data KEDA of the mechanism was carried 
out. The results obtained are represented graphically as shown in 
Figures 5 to 16. 

Figure 5 represents the fluctuation of angular velocity 
© of the crank. The desired average speed of the crank is 21.0 
rad/ sec while the actual average speed of the crank is allowed to 
fall in the region represented in between the lines 1 and 2 (FIG-.5)? 
permitting an error of 0.5^. 

Figure 6 represents the variation of slider displacement 
X and "rigid" displacement X,, of coupler - end C with respect to 

U 

crank rotation 0. In the absence of spring and with rigid links, 
curve 2 would have represented the variation of slider displacement 
X. Therefore, the vertical ordinate between curves 1 aid 2 at any 
value of 0 represents the combined effect of inclusion of the 
spring and inherent elasticity of links of the mechanism on slider 
displacement X at that position of the crank. It may be noted 
that the characteristics of the mechanism are calculated taking 
into account the variation of crank speed © due to the changing 
inertia of the system with crank rotation. 
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Figure 7 represents variation of angular velocity of the 

coupler 0 , vhile Figure 8 represents the variation of slider velo- 

* 

city Z and x and y - directional velocities of the centre of gra- 
vity of the coupler with respect to 6. As nay be expected, the 
slider velocity is not zero at crank angles 0° , 180° and 360°. 

Figure 9 represents the variation of angular acceleration 

ox crank aid coupler with respect to 0. Refering to Figure 5, 

• • 

since the crank - speed decreases from © = 0° to © = 48°, © is 

• • 

negative in this range. Between © = 48° and © = 130°, © is posi- 

tive as © continuously increases in this range. Similarly between 

© = 130° and © = 275°, © is negative aid from © = 275° to © = 340° 

* # • • 

9 is positive. From 340° to 360° of crank rotation, © is again 

negative. 

Variations of linear accelerations of centre of gravity 
of the coupler in x aid y-directions and slider acceleration with 
© are shown in Figure 10. 

The variation of elastic deformation X at coupler end 

e 

G with © is represented in Figure 11, while Figure 12 aid 13 re- 
present how the remaining element deflections vaiy with ©. The 

elastic deformation X (FIG. 1l) when added to 'rigid' displace- 
— e 

ment X (FIG. 6) will yield the actual displacement X of the 
coupler end G. The variation of the length of spring (X - X ) 
with respect to © is shown in Figure 16, wiile the variations of 
inertia forces with © are plotted in Figure 14 and 15. 


, , -25631 * 

If*', 

ICi 

8n. 
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Befering to Mgures 9, 10, 11, 14, 15 end 16 it may be 

noticed that there is a change in the nature of the curves of 

• * • * 

9, X, X , P , P„ and (X - X ) in between © = 160° aid © = 260°. 
e J *+ o 

mis may be attributed to the nature of variation of X_, aid 
, . e 

hence (X - X^) with respect to ©. The clastic deformation X^ 

e 

varies between 160° aid 260° of crank rotation in such a way that 

(X - Xq) is affected in this range of 0 as shown in Figure 16. 

The variation of (X - X ) directly affects inertia force P. 

u 4 

(equation 2.30) and slider acceleration X (equation 2.l), which 

♦* 

in turn affect the input crank acceleration 0 and inertia torque 



The variations of slider displacement, velocity and 
acceleration with input angle © as obtained by Davidson ^5 J for a 
schematically identical mechanism, are shown in Figure 17. Compa- 
ring his results with those obtained in the present analysis it 
was found that as may bo expected, they do not tally. The diffe- 
rence in the results of slider acceleration is much more pronounced 
in compaxision to the differences in displacement and velocity 
results. As per Davidson's results the slider moves 4.0 in. as 
the crank rotates from 180° to 300° while the present analysis 
gives this distance as 3.99 in. Davidson predicts the slider velo- 
city at © = 180° as 0.0 in./sec. , at © = 240° as 42,3 in/sec aid 
at © = 300° as 66.0 in/sec. The corresponding values obtained by 
the present analysis are 2,269, 39.352 and 67.8 in/sec respectively. 
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Since Davidson did not consider the effect of inherent 
elasticity s the fluctuation in driving link speed and the mass of 
crank and coupler (which have been accounted in the present anal y- 

sis) the present results for Z are much different from one obtained 

#• 

by him. He obtains values of Z at ISO 0 , 250° aid 300° of crank 
rotation as 1900, 0.0 and - 900 in/sec^ approximately as against 

O 

964.4, 740.0 and - 300.0 in/sec obtained in the present analysis, 
3.2 Examples 2 and 3 

The mechaiisms studied in exaaples 2 and 3 are schema- 
tically identical to those synthesized by Liniecki : 26~j . In 
Liniecki* s first model the radius of crank, length of coupler and 
the amount of offset are taken to be 5.889 in., 19.164 in. and 
0.688 in. vfaile in the second model those values are taken as 
5.933 in., 26.961 in. and 2.59 in. respectively. In both the models 
following values have been adopted by Liniecki : 

-1 2 

Mass of the slider = 0.52 Ibf in see 

Moment of inertia 2 

of crank about its = 24.0 lbf in sec 

axis of rotation 

Driving Torque 

Constant =1.8 

Average speed of 

crank _iX = 26.2 rad/sec 

Since Liniecki has considered rigidly attached slider, 
the stiffness of the spring in the present case is taken to be 
very high (52000.0 lbf/in.). Liniecki has not considered mass of 
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the coupler separately (he considers equivalent moment ox inertia 
of rotating masses and equivalent reciprocating muss in the mecha- 
nism) and has assumed '’rigid” links. The following additional 
infcnnaticn is assumed for both tile models examined here, 

2 

Area of cross-section of crank = 2.0 in. 

2 

.Area of cross-section of coupler = 1*5 in. 

6 2 

Modulus of Elasticity E = 30 * 10 lbf/in. 

Moment of inertia of coupler 2 

about the transverse axis = 5.1 lbf in. sec. 

passing through its centre of 
gravity. 

The mass of the coupler is assumed to be half that of the slider. 

3,2.1 Results and Discussion of Bxaaple Problems 2 and 3 

The KEDA of the mechanisms of examples 2 and 3 was per- 
formed, The difference between the actual velocity of slider end 
the desired velocity [j26j w as obtained for every 20° of crank 
rotation. Finally this quantity was plotted in Figure 18 against 
crank rotation 0. 

Liniecki synthesizes the first mechanism by making the 
error between the actual and the desired velocities of the slider 
zero at 60° , 80° and 120° of crank rotation. The second mechanism 
is synthesized by making this error zero at 60°, 160° and 240° of 
crank rotation. But the results of KEDA of both the mechanisms 
predict that the velocities at these precision points will not bet 
zero (Fig* 18) if elasticity is also taken into account. This 
error at precision points varies from 3% ^o 10/&. 





CHAPTER IV 


CONCLUSIONS AND SCOPE FOB WETHER TORE 

Tiie importance of KEDA of mechanisms has long been 
recognized by the researchers. As it is clear from the exanple 
problems worked out, this analysis is essential for the mecha- 
nisms where accuracy is main objective, since the elastic 
deformations may cause appreciable error in mechanism charac- 
teristics. 

However, in the present dissertation the effects of 
clearances and tolerances, danping characteristics and friction 
in joints have not been taken into account. An analysis con- 
sidering all these effects together would be an ideal one and 
will predict the actual characteristics of the mechanism much 
more closely. Efforts in fixture may be directed to formulate 
a mathematical model of mechanisms considering all these effects 
and to device suitable solution techniques. 
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132 FORMAT ( 50X, *VALUE OF XPR IM*//6 ( i r ‘X , F* 5 . 3, 5 ( 5/ , F ±5. .? ) / ) 3 

DO 200 1=1, 16i 

23 J XDOT ( I ) =Z Z ( 1, 3 ) *S’3RT ( ZZ ! I ,2)3 
PR I NT 1 03, ( XDClT ( I ) , I = J » 131 * 3, T } 

*03 FORMAT! 5DX, *VALUE Jr XD0T*//6 ( 10.. » Fj. 5. 3, 5 ( 5X , Fi5 . i ) / ) ) 
PRINT256 


t 




• . I j r I . . + 3 


256 FORMAT ( 5'JXj *C.J ! , ’ 1 / c r'3 ' r- 
83 3 CONTI NUi- 

1 = 

X = * ' »■ 

P3=R**2/( 4. C * 4 L ) 

P4=A**2/( 2. U-*^L 3 
P7=R*A/AL 
P1>R/AL 
3 1 = 1 + .. 

I M I TA ( I ) = o 
11=1+1 
1 J = I +2 
xX=2.0*X 

ANGVl'U I)=SgR1 (V(I) } 

P5=WN**2/ { Wi i**i-Y ( I } ) 

P6=WN**2/ (WM c *2-4 . * Y { I }) 

P8=AL* ( 1. 0- { A/ Al + R* 3 I j{ •')/,- L 3**1 j.Z 

P9=P8**2 

Pi *= ( A + R*SIM X ) ) *(R*COu(X ) )**•' 

P12 = R*SIM (X 3 + • . 5 * P .. .*{ . 5 ~R*$Iw( >() + ’.*: JS < X) 5 
P13=R*C0S (X 3+0 .5* PI , *{ H^C-.'S ( XV 3 - A* _ i ;[/.)) 
Z=aRSIN(A/AL+R*SI 4( s ;)/.\l) 

P14=F*A0 

Pi5=h*i; 

Pi6=2*4G 

P17 = SIN(X) + T^CZ)*C JSC i) 
Pi8=C0SU)-T*,j(Z)*S! - .(,:} 

P19=C0S (X )-.5*TANCZ 3*81 ,( O 
P20=S I N { X )+ .5*TaN (Z )*C0:>( ' 3 
P2 *. = T AN ( Z ) / CCS { Z ) 

P22-1. ?/Ci'jS (Z)**2 

P?3=Pi3**I 

P24= Pi 7*4=2 

3 Ri = ; . 3 *AL**^*P* 7 /P. T 
BB2i=-S IN ( X 3*P_S/CP (Z3 3 

B822=AL*Al *CU$(X)*P„7/{ .*P: 5*C r iSCZ) ) 
BB^3=P21/ PI b 
3 B 2 = S R 2 1 + <3& ? 2 +B B 2 

8B3i=AL*P19*P„j/F.,4 + .5* >,L *P 6 

BB32=AL*AL*AL*Pi7*P :n/{. .A*? *3 3 
BB 3 =- B8 3 1 -R B 3 2 

8B4=AL*P23/P*‘r+AL*P2'‘/P .6+AL **^* P? '>/ { 3 . *Pi53 
BB51=AL*Pifi*SIN(X)/!2.* D . \) 

8B52=AL*A L* AL*P17*C'J3 C X 3/ ( 6. *Pj.5 3 

B85=8352-BB5i 

XC ( I 3=R*C0S (X3+P8 

ALP l=AL/( E*AO) 

ALP2=AL**3/ (i.O*E*I _) 

ALP3=AL**2/ (2.0*£*I . 3 
ALP4=AL/C 3*113 
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4L P? = AL*P 1L / ( . * ; 3 ) 

ALP 6 =Sir;{ Z)*p_„/( * ,) 

ALP7 = --»L*P? . ; ;• 

ALP- = 0 L/{ : ? . 3* ■.« I ✓ } + 7 ; {z )-<■-= /{:]_• -?■ . .) 

AL P9=-0 .5 -ALP... 

ALP., - A L * { P o S ( Z ) H * . ' / -t a * I + ~ i * v i I = l L ) :: 3 / ( - ■ • * * I 5 

*lpj : =al**s*si;i< z. *zj/d.. • . *z )/{•,. *z*'G > 

alp:2=4l*si ,;( z ) +*_/ ; ' . „* -=^ 7 >+, l*~-' *c kj{ z . j*=-*i-) 

BL=— SI<%! (X )/{ j&L-^CCCC ZS ) 

82=-P?.9 
fc3 = P^ 

B<-=-o.5*si . { ; 
b 5 = C 0 S ( X ) /{ U*CuS(Z S ) 

B 6 =-P'Z- 
B7=P i 7 

68 = 0. 5*C3 3( K) 

8 B 6 =A LP 3 * Zb 

dB7=ALP 3 * B 1 ** • + ^L r '--r 7 *.: _ + i_ r. ' 

B B 6 = A L P 1 * 8 1 * B _ + * L " 7 * . 5 * » < j + -• L p > 

BB 9=AL P 1*D6 #L i + A L P Z ! ’ ! * r ^ ~ + A L P c 
861 .=ALP1*C_*L,'V+ALP_*B.;*‘' ' 

8 B x 1 = ALP3 *S o 

B812=ALPl*3i*3Z+ALP 2*3 5 *b 6+- 1 = 3 
BB3 3=ALPl*Bi**2+ALP , + ALP . 

BB 3.4= A L PI *5 2* 33 + ALP '.* 86 * 8 ' +ALP? 

BBx5 = ALPl*5<i* -V + A L P ! s - : V c v L . •' + 4 L P - . 

BB 16 = AL P3 *3 C 

SB 17 = A L PI *b i *34 + ALP;; *>35 ** 

BBx 8 =ALPl *B i*R 4 + ALP ■ ’ * d 6 * n ,> + A LP* 

BBi9 = ALPl*3'* r 34 + ALP P5?*-., 

BB 20 = AL »i *B '•: * * 2 + A L D 2* . 3 + *_ +A L P . 3 
IFURE.GE.i )GU Ti! 962 
XCF( I )=XC ( I ) 

962 CONTINUE 

XF DOT ( I ) = - A NijVlL ( 1 ) *(R*pv*SI . f X } +? *P 5*C r JS { X )+2 .*Pb*P6*S 1 M XX ) ) 
ACCXG ( I ) = -P.2*ALPMM I )-P..;*Y{I ) 

ACCYG ( I )= .5*R*(C0S ( / ) * ' L PHA ( I 1 - SI .( A)*Y( I) ) 

BETA{ I )=(R*CUS(X)*AL=>h; (1 ) + (i J * / P9-R*S IN ( X ) )*Y( I) )/Pq 
Q1 ( I )=0I*ALPH4( I ) 

Q2 ( I ) =GI* 8 ZTA{ I ) 

Q3 { I )=0M2*ACCXG( I ) 

04(1 )=S*( ZZ (1,1 )-XCF( I ) ) 

Q5 ( I ) =0M2*ACCYG( I ) 

XCD 6 L ( I ) =BBi*Ul ( I ) +882*91 ( I ) +BC' *•» A I) +BB4*C4 { I )+RB 3*w5 ( I ) 

XCF ( I )=XC(I ) +XCDE L( I) 

555 CONTINUE 

IMITACI )=IMITA( I)+l 

PHYDEL ( I ) =B 86 *C„ ( I ) +937*0 . ( I ) +- n 2* ( I )+B89*Q4 ( 1 5 +BBiV*Q5 ( i ) 
XGDE L { I ) = SB I ' *0" ( I ) +RBJ '*22' I ) + •' cl *03 ( I )+BBi4*Q4( I )+3Bi3*Q5( I > 
YGDEL ( I ) =83 1 6 *G 1(1) +13„7*02 ( I)+oB. *Q3 ( 1 ) +5 5 ‘.9*2- (13 +3329*3 5( 1 } 




IF ( I . EQ.M ) SC f >"!’ 5 

HHi={PHYO=L (ii )-PH Y - L( (3 3/T 

HH 2= ( X 3DE L ( x I } — •' >3 Z L f I ) )/' 

HHJ=(Y9PEL( II )-YG.:i {!))/- 
IF ( I.FO.UM I J = ; 

HHs-= ( PHYDFL ( I J )-2 ; (II3+PFY :LtI ))/~t*2 

HH 5= ( XGDF L ( i J ) - 2 . i * - _, L , .>_( < | } + vr, , t j_ f j j ) /j**; 

HH6=(YGDEL( *YSD l ( I I ) +Y^": tL { I ) 5 /T**2 

08 ~T A ( I ) = HH4*Y ( t )+Hr, * £ _ P L ( I } 

DACCXG { I ) =H ! ;3 v { J ) + , 1-2* . _ ~ k • I ) 

DACC YG( I ) =Hh*-, * ( ( I ) + ,l -’HA { t ) 

GO TO 552 
551 CONTINUE 

D3 ET A ( I ) = 03 ~T ,( ; 5 
DACCXG { I) =D~XC>;G( „} 

DACCYGt I)=DACCYG(.. ) 

55 2 CD NT I Nllfc 

BETA { I ) -D8ETA (I ) + (.U..ns(, ) *fi LPH ft (I i+i^ 3/ ?'*-?.#$ I I { X ) ) *Y ( 1) ) / = , 

ACCXG ( I ) = -PU'*ALPf 4(13- >_\*Y( I )▼£•’, CCXGt I 5 

ACCYG ( I ) =~DACCYG{ I 3+ ( COCO } *A_PHA (13-SI j(X)*Y( I) ) 

OKI ) =GI * AL PH~ { I 3 
Q2( I ) =GI*BET A ( I ) 

U3 (I ) =0M2*ACC*G( I 3 
04(1 )=S*( Z2 (1,3 )-)XFC 1 3 ) 

Q5 ( I )=0M2*ACCYG( I 3 

XC DEL ( I } = P.R ..*CJ3 ( I ) + O r 2- l( I H-Pf.li* ■). ( I )+& c ,4*04(I > + 535*25(1} 

XC FI ( I )=XC( I ) + /XC; L ( I ) 

HHK=(XCri(I l-YCFU) ) / XC F ( ; ) 

XCF{ I ) =XC Fa. ( I } 

I F ( ABS ( HHK ) . LT . - PS-+ ) GO 12 sS*. 

GO TO 555 
554 CONTINUE 

IF ( I - EQ.M >3fJ TO = . 

DPHDOT ( I)=HH„*AnGV2L(I 3 
OXGDOT( I) =HH2*4nGV2L ( I 3 
DYGDDT ( I ) =HH3*ANGV-L( I 3 
GO TO 5 57 

556 CONTINUE 

DPHDOT ( I ) =DPHDOT ( . ) 

DXGDOT ( I ) =DXGD(jT ( i ) 

DYGDOT ( I) =D YGDOT ( 2 ) 

557 CONTINUE 

PH YDOT ( I) PHDOT { I ) + Pj.**C’3S ( X } *AAG / L L ( I J/CUS C Zl 
YGDOT ( I ) =-DYGDUT l I 3+ .5*R*CUS ( X 3 * A J D Vl L ( I 3 

XGDUT ( I 3 =i)XGDU T ( I )- A tGVCL ( I ) *( R* SI . ( X 3 +0 . 5*P.l* (0. 5*R*SIN( X X) +A* 
1C0S(X) ) ) 

X= X+ T 

IF(I-LT-M) GO TU 1 

PRINT638, ( I N I T A { I), I = it ,4 .10) 

688 FORHAT(10X,*NO OF IT' RATI- NS BuFCR. CDNV GGP.C: */!X , 2G I 6 3 
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PR I NT 201 , (X D j'JFL ( 15 , 1 = . • , r, \ 

201 FOR lAT(5DX, *V-L'j : ,P 

PRI-iT? 2 » (4 .G'J „L ( i ) , = t , 

232 FQR 4 AT (/5 ’ J : »‘ L 

PRlNTPOi, ucc '(ii,J= 

^■•0-' FORMA T ( 50 A, +SL1 ’Jl R C>V/-( . . , f 

PRI <T2.:4, (XFXiT( I) , 7 = i , _ 3 

234 F3RM^T(/5 X,*wMLL„ .'•* ,,'j “ .;<// < 

PRI NT 20 5 1 (£LPH*.{ !),:=, ", ) 

205 FORMAT ( /50X » * A «r} »C'„*//_(. ' , ’ 5, 

P R I it T 2 6 » ( B - T A { I } , I = ) 

236 F0RriAT(/5 ■*, !‘V , ,L , 2 ; - t *//-,_ 

PRI FIT 2.07* ( 6 CC \0 ( I . * _ : > 


, 5 f X f P ) / ) 3 
rf ’ - { . , 5 ( ’ X , F .. 3 . 0 3 / ) 3 

. » = ( 5 v , F . « . 3 ) / ) ) 

- r • ' ) " ( jAfF, j , l ) / } ) 

r ! / 3 5 

, > X , F 3/3 ) 


5 ' 5 Zt / r f* 


237 

FORMAT t/5 OX, *V/LL. - 


: 'i ■'// - 

( , 

■ , :. ( 

s A , F . 


PRINT’ >73, ( ACCYJ( I } , 

i = - f 

• - . ? 3 




37 > 

FORMAT ( /5 ' X , -r i,C C YG 


■’-// A 

A f ' 

. . , = (5X 

, F_ 5 . 


PRI t'.T 212, (0. ( i 3 , I = .. , 

* t 

. 3 




212 

FORMAT (/5 OX, *VmLI. 

'It 

0 

, # C 


X , F 1 h 


PR I NT 213, !wx( i 3 , I= 4 , 

'■J . 9 

. ) 




213 

FORMAT ( / 5 'X , *\/„tL - 

• IF i 

’ _*//5{ 


3 . : , 5 ( f 

.:,f 1? 


PRI NT 214, (C. ( I 3,1 = 1, 

-- i 

3) 




214 

FORMAT { /p X,*VML. 

jF ■ 

■ #//t.{ 


: . 3 , 3 S c 

: , f i ' 


PRlNT?i5, (Q4( I 3, I=> ) 

215 FORMAT ( /5 ,Xt*VALL , p - •' F 7 *//£,< 
PRINT2152,(*j5(I),; = . , 1 T , 03 
: 15 2 FORMAT ( 50X, *VaLGE . IF ,;v// 0 {. 


. : 3 / 3 3 
. j ) / 3 5 
. i) / 3 3 
5.I,5(5X»FI5.o)/) 3 

. , J ( = : » f 5 . i ) / ) ) 


PRINT716, (XCO ;i( I), 1 = 1, V.,: 3 

216 FORMAT ( /50X ,* VAUJ2 -:F 'O-uL -- / /b { ,_0X , F! 5 . 2 , 5 ( s v f F * 5 . 0 } / 3 ) 

PR i NT 217, { X C ( I 3,1 = 1,*:..,. '■) 

217 FORMAK /5 X,*VALL e • .F SC*ffi { ”, x , F _ F . o , 5 { 5Y , F 15 . 3 3 / 3 3 
PR I NT 2 19, (XCF{ I ),I=. ) 

219 FORMAT { /5 OX , *F I ivlAL >'CF*//6( ..O'. , F. 3..- , a ( 5X , F_ 5 . 8 ) / ) 3 
PRINT220, (PHYDiLU 3 ,I=. ,1 _, ; A ) 

22 w F0RMAT(5\,X,*PHYCSL*//M. . , 5 ( , F 1 5 . ? 3 / 3 3 


PRI NT 221, (XGO-Lt I), I = A , ... . ,1 3 

221 F0RMAT(5DX,*XbD>-L*//6(.cy,F<.5.<.,5( .,.= -5.2}/)) 
PR I NT 22 2, (YGU1L( I 3 , I = 1 , „ o. , 1 ., ) 

222 FORMAT(50X,*YGD:-L*//6(* FI 5. « , 5{ • F. 5.G) /) ) 
PRINT223, ( PHYDGT ( I ),i=_,_^,_ 3 

223 FURMAT{50X,*PHYOCT*//b{ 9 Y » F . a . " , a ( 5 '• , F j. 5 . <3 ) / 3 ) 
PR I MT224, { XGDOT ( I ) , 1 = 1 , _ 0 „ , „ 0 ) 

224 F0RMAT(53X,*XGD0T*//M1 < , FI b. ft , 5( ~ X, F . 5 . 8 3 / ) ) 
PRI NT 225, (Y GD'JT (13,1 = 1,--. ,1 3 

225 FORMAT(50X,*YGD:iT*//M„0X,FI,5.:.,E{ .,F*5.«)/)> 

PRINT226, (DPHDOT(I) ,I = . ,_ 31, , 3 

226 FORMAT ( 50X, *DPHCCT*//6( . : : , F .5 . 3 ,5 ( 5X , F 15 .8 3 / }) 

PRINT227, (DXGDUT(I) ,1 = 1, . 3 

227 FORMAT ( 53 X, *DXGCCT* //6 ( . JX, F 15 . 3 ,5 { 5Y , F15 . A 3 / ) ) 
PRINT228, ( D YGOGT ( 1 3 , 1=.,--. ., 3. 3 

228 F0RMAT(5DX,*DYGDCT*//6(„ v , F 15. t , 5 S 5X , F 15 . 8 3 / ) ) 




PRINT229, { 3 hr~ A ( I 5,1=1, . , . , j 

229 FORMAT ( 52X, *DC iTA-*/ /*, ( . " , F_ r' . « , 5 ( ' / , F . 5 . ‘1 ) / ) ! 

PR1NT2S'*, (D«CCXG{ I ) , l = .,.v', } 

230 FORMAT ( 50 X, *DACCXV*//6 { ; 'r. 1 , - t F'5. )/}} 

PRINT231, (DACCYS(I) , J=, , , ) 

23 i FORMAT { 50 X, *DACCYG*//f { . , F _ S . ;»"* 3 •" . » F _ 5 • ’-,)/) } 
DO o 3 7 1 = 1* ju ^ „ 

HK = (ZZU,1)-<P(I))/2Z(I,.) 

I F { A b S ( HK ) . G _ . 0 . 1 ) GO ~j 

837 CONTINUE 
PR I NT 339 

839 FORMAT { 50 X,*PRG3L 1 HAG CO, .7 - AG j, , H* 
Re TURN 

83 8 DO 84.* I = ; , 3 . 

84" XD { I > =ZZ ( I,; ) 

IRc=IRc.+i 
PRINT961, i R f 

9Si FORM/ T ( /5 X,*R.F11 Q = *,.'/* \ , * : , ..n* } // ) 

MM= , 

836 CONTI -!iif 
MM = Mfi + j» 

K= 2 
X = ® ■ 

DO 809 1=1, a 

809 ZX (I)=ZZ( :<* 1 } 

DO Gi.' J K = a , N 
KK=JK+i 
KK J = KK+i 
X], =X+ : . 3*T 
X 2 = X + T 
XCFF=XCF( JK) 

XCFP= (XCDcL (KK)-XCD L(JK) )/~ 

ZS ( 1 ) =PHYDOT ( JK 5 
Z S ( 2 ) =XGDOT ( JK ) 

ZS ( 3 ) =YGDOT ( JK ) 

1 ZS(4)=3FTA( JK) 

ZS ( 5 ) = ACC XG ( JK 5 
ZS { 6 ) =ACC YG ( JK) 

CALL FUNCt X, ZX»DZX ,K,XCFF) 

ALPHA ( JK)=DZX(2)/2.'J 
XD OUB L ( J K ) = D Z X ( 3 ) 

Y ( JK ) =ZX ( 2 ) 

ACCX ( JK ) = ZX ( i ) *ALPHf { JK ) + Z; { 2 ) # ,/ CO •/■8L ( JK } 

DO 82" 1 = 1,3 

DZY1 ( I ) =D ZX ( I ) 

820 ZY( 1 )=ZX( I) +’ .5*T*DZ<{ I ) 

XCFF=0.5*(XCF( JK)+XC p (KK) ) 

XCFP="*.5*(XCFP+{ XCD ;L ( KKJ J-XCuaL (KO)/T) 

ZS«1) = j. 5*( PHYDOT { J K ) + a HYDuT ( KK ) ) 

ZS ( 2) =w .5*( XGDOT ( JK ) +XGPO T {KK ) ) 
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ZS (3 ) = , ,5*( VG3L T ( JK } + Y'j ._T { <.< } ) 

L S(< } = '.5*(H"-H(JK)+i r T ) 

ZS ( 3 ) =0 . 5 *( ACGX j ( Jr } + _ ' - - / / - ) i 
ZS (*.) = .'.?*( - C C Y 3 { •<) + {^ 5 ) 

CALL FUNC ( X. ,ZY»~Z, , < , . "Fj 
DU L ‘ 3 I = „ , . 

DZ Y2 ( I ) =3 ZX ( I ) 

82 5 ZY ( I ) =ZX< I)+ .p*T**, 4 (I) 

XC FF= .5*(XCF(JK)+> W F ( .</. } J 
CALL FljMC ( X „ , Z Y » 1Z ,> , J , i - z F) 

DO P.30 ! = *„,?. 

DZ Y~ ( I } =3 ZX { I ) 

83 - ZY ( I )=ZX( I) +T*:ZX( I ) 

XCFF = XCF( KK ) 

XCFP=(XC3LL(KKJ)-)C.-.L(. ) )/ ; 

ZS (I ) =PHYD3T( KK ) 

ZS (2)=XGD0T ( K K ) 

ZS ( •> ) =YG30T ( KK ) 

ZS (4) =3cT a{ KK) 

ZS ( 5 ) =4CCx3 ! KK ) 

ZS ( 6 J =ACC Y3 { KK ) 

CALL FUMCf >Z, ZY S JZ ',K, C F F ) 

DO 83 5 ! = 

DZY4( I ) =DZX { 1 ) 

ZX { I ) =ZX! I) +T*(DZY: ( I ) + ■ . 43 ZY • c !) + ''. *DZY • f I ) + DZYt( i )) / 6.0 
835 ZZ(KK,I) = ZXU) 

X = X + T 

81 j CONTINUE 
MK=3 
NK=? 

GO TO 824 
END 
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*IBFTu FUNC 

SUBROUT IM f FU'jC ( \ T Z i OZ ? F , aC fp J 

* * * * * * * ** $ * * * * * * ^ ^ ,; c 4 , 4 # , ;s * v v ^ ^ i: 

THIS SURROJTI ;= 3V aLM^T.-S H- F 'f,C I 
APPLYING RU jG iT 7 . iThF.'^F ’ T 

COMMON/TAJ/ZS 
COMMON/ SM IT A/XCFP 

COM^QN/khARl/ r» AL , A, ;• I } ,1 , *k', ,,,v 

DIMtNSlOM ZSC6) 

DI MON SION Z {.-,)} C ? ( - ) 

XX=2 • *X 

P 1 = ( A+R^S IM (X ) )/£L 
P2=A/AL 
P3 = R/AL 


7 A L . > K ; S 1 1 .1 _ ^ 
S “ L- ’ 7 I M .'F A 


- - t - 1 - - 


P4 = P J. *#2 
P5=i . u-P4 
P6=si n ( x) **:: 
p? =cos ( x) **r 
p s=cns«x) ** j 

P9=L,I+GI*Pb**Z*P7/P - 

P I 0 = 0K>2*R 2 * ( P 6+ j • I f * D i 1 -■ 5 :Ji r+ ;■ . -i : C I { ’ / ) *P 1 ) 

SS**P9+P1. 

pi i =gi*P3**2* ( 2 . ■:*? yrpj.-.p -si ■ ( ■ x ;• ? 3 ) /? 5**2 
P12=OM2*R**2 : M SIM \X) + . 5*p~ *0 - *p_- j. 25*5 I I ( >' A } *P ) 
Pi 3=OM2*R **«!*{ COS ( X <) *P - + 31 , { \ X) *CHS(X ) ) 

SS2=P1 1 +P I2+P • 3 
SS3=-R*C0S(X)-AL*?3** . : 

SS4=R*S I N (X }+R*CCS( a)*p. 

SS5 = WN**2*{ R*LGS If ) + ■ L ) 

H5 ?’=OM2 *R ** Z / 4- , . 

H5*=H54*C0S <X )**2 


H52=H5 i*S IM (XX) 
SSju=SSi+H5i 
SS2 = SS2-H 5 ? 


DZ (1 )=Z(3 ) 

P14=A1-S*(Z{:)+$S')*(Z( J+SS« )-Z (T >*{ A2+Q.5*$S2) 
UZ (2)=2.0*( P14-Lfi*Z(M*( S S 3 - Z { 1 I/SSi 

DZ(3)=(SS5-Z(i)*W\**.- . 5-Z C 3 ) -7 £ [ _))/Z(2) 

IF { K.FQ.l ) RETURN 
SS4=SS4-XCFP 
SS3=-XCFF 
SS5=W-N**2*XCFF 

P14=A1-S*(Z (l)+SS,f)*(Z( )+SS -O-Z { ? 5 2 

P15=GI*ZS (i )*ZS(4) + u.-2*(ZS(l)*Zi (5 ) + ZS ( ':• ) *ZS ( 6 ) ) 
P16=-Pi 5/ SQRT ( Z { 2 ) ) 

P17=OMl*Z(3)*(Z(l)*W 3S5 } 

DZ (2)=2.D*( PI 4 + P 1 6+ P \1 ) /HI 

DZ { H ) = ( SS 5- Z ( 1 )*WN**i-0.5*Z( 3 ) Z { _ ) ) /Z { 2 ) 

RETURN 

END 


•i NIL; 
iLY.ii 





*tNTRY 


'i' -J* -V “J * 't' *1' -<* -T' 

- : p ! i t :a‘T\ * 

4- ** ' ' *-< *• J. b V aji«r •**. V J*. .i. 

"V” -v '«*• 'T" *•** -i e- *r ** *r 'c- *f- '!* -»•• 


12.55 

445 . 'J 

' ® 

1.0 

: . 1 53 

. * l , 

v5 • 3 i5 j 

. o46 9 

0.0 

32.78 

700.0 

-0.10 

24. '5 

5.1 

* > 

v . 3 3 3 9 

' .188 

* * 

24.9 

7 ’ v.O 

-0*5 

24.0 

5.1 

0 . C "V 

0.3334 

• .188 

„ • „* 



4.38 
>-.5 
o . u 

V . 9 VI 


_ 9 • ± h -* 


x • , 

U.'JUb 


0 . 6 b J 
2 * - 
■ * j.i. 


1\> IV 
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